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We examine the entanglement between two qubits, supposed to be remotely 
^ 1 located and driven by independent quantized optical fields. No interaction 

is allowed between the qubits, but their degree of entanglement changes as 



a function of time. We report a collapse and revival of entanglement that is 
^ | similar to the collapse and revival of single-atom properties in cavity QED. 
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^ . Control of the evolution of qubit entanglement and the time-dependent behavior of qubit 

OO ■ 

q ■ pairs in networks is relevant for quantum computing and cryptography. Entanglement will 

t> 
o 



be stored in quantum memory registers (see conceptual sketch in Fig. 1) for eventual use in 
some form of quantum communication. The qubits must be controlled in some way externally, 
and we are interested here in the response of the entanglement of a pair of stored qubits to 
quantized optical control fields. 

What we can call "pure storage" requires mutual isolation and non-interaction between 
qubits. Two-qubit evolution has been studied without qubit isolation and usually allowing 
or relying on mutual interactions to produce entanglement dynamics. Kim, et al., showed [1] 
that an incoherent thermal field can create entanglement between two such qubits. Entan- 
glement transfer between two qubits and two separate fields was examined by Zhou and 
Wang [2]. In their treatment the quantum field was weak rather than strong. The evolution 
of entanglement in a qubit-field system, where the qubit and the field start from mixed states 
was examined by Rendell and Rajagopal [3]. They aimed to calculate the entanglement em- 
bedded in the full system, and because of the lack of an entanglement measure for 2 x oo 
systems they calculated a lower bound for the concurrence instead. 
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Fig. 1. Sketch indicating non-interacting qubits in a quantum storage network. Dashed line 
indicates two are entangled. 

Other studies [4, 5] have shown that two isolated qubits can exhibit periodic fluctuations 
in their entanglement in the form of early-stage decoherence (ESD - also referred to as 
entanglement sudden death) [6] when the qubits are modelled as "controlled" locally by 
interaction with only single photons. However, manageable control fields are better modelled 
as containing many photons. Here we retain a quantum picture of two many-photon well- 
phased control fields by using a coherent state description of them with a large mean photon 
number n ^> 1. 

For our calculations we take single-mode control fields. Each field is assumed, for simplicity, 
to be exactly resonant with the flip transition between the ground \g) and excited |e) states 
of the qubit that it addresses. The well-known Jaynes-Cummings (JC) interaction [7] is then 
relevant at each qubit site (labelled i = 1,2), and the interaction Hamiltonian is given by: 

H l = n 9(ai°t + a l°i), (!) 

i=l,2 

where a, and a\ are the photon annihilation and creation operators for site i, and af and 
o~[ are the raising and lowering Pauli matrices for atom i, and g is the coupling constant 
between atoms and fields, taken the same for both for greatest simplicity hereafter. For the 
two coherent state fields we take the same n for simplicity, and we assume initial entanglement 
in the form of a familiar Bell State: 

|*(0)> = (\eg) + \ge))/V2. (2) 

Since JC time evolution is unitary, no information can be truly lost in evolution via Hi, 
but the infinite range of photon numbers in a coherent state brings aspects of open-system 
theory into play. However, the quasi-classical nature of practically available control fields 
suggests that we not expect their quantum characteristics to be dominant. Thus we will 
trace out the fields and follow only the qubit entanglements. 

In this discussion we will use Wootters' concurrence [9] as our entanglement measure, 
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which is given by 

C(p) = max{0, \f\i - \J~\ 2 - yAjj - ^/V}-, 
where the quantities \ are the eigenvalues in decreasing order of the matrix 



C = p{<7y®Py)p*{<7y 



(3) 



(4) 



Here p is our two-qubit reduced-state density matrix, p* denotes the complex conjugation of 
p in the standard basis, and a y is the Pauli matrix expressed in the same basis. 

The photon number in a coherent field is Poisson distributed and relatively tightly centered 
around n when n 3> 1. This suggests a shortcut approximation, to be checked numerically, 
in which we represent the field density matrix as a Fock state having photon number equal 
to n. An important simplification occurs in taking the same n for both control fields. Then 
the initial field state is \n) ® \n) and the reduced density matrix for the qubits is: 
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where we have used the standard two-qubit basis [ee, eg, ge, gg\. The elements indicated 
by x are zero because of the equal-n simplification. Thus, under the assumptions mentioned, 
p is of X-type (see [8]). 

For an X-type p, Eq.(3) turns into: 



C(p) = 2 max[ 0, 



ad 



(6) 



The control fields induce growth in time of the elements a, d, which are the only ones not 
already present in the original maximally entangled state (2). Their growth and any decline 
of z will cause entanglement to decrease. 

Having used the Fock state shortcut to obtain (5), we avoid using it further now and 
calculate the elements z, a, d for the coherent state. We introduce the Poisson number 



distribution by the coherent-state amplitude measure A n = e 
we obtain a doubly infinite series summation: 



-\°\ 2 /2 n n 



a n / Vn\, where |a| 



n, 



n,m 

+ A n A n _2A m A rn+ 2S n S n -.iS rn+ iS rn+ 2 
— A n A n ^iA m A m j r iS n C n ^iS m j r iC m j r 2^, 

where C n = cos(gt^/n) and S n = sm(gty/n). 



(7) 
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Similarly the series summations for a and d are; 



+ AnAn+iArnAm-iSn+iCn+iSmCrn 

+ A n A n _iyl m A m _|_iS' ri C , n S' m+ iC , m _|_i| 



and 



+ ^n^n+l^m^m-l'S'n+lCn+l'S'mCm 

+ ^n^n-l^m^m+l5'nCn'S'm+lCm+l}- ^ 

The infinite extent of these summations of course reflects the fact that we have coupled 
the qubits to an infinite state space. The sums cannot be analytically completed, but in our 
calculations we choose a = 10, i.e., n = 100 ^> 1, so we can obtain good approximations if 
we use Stirling's formula for n!, 

n\ = v / 2^n n e" n , (10) 

and Euler's formula to approximate the terms in the summations above by integrals. If we 
approximate the terms like A n A n+ iA m A m _i with A 2 n A 2 w which introduces an error of order 
1/n near the Poisson peaks n m m « n, helpful cancellations can be identified, and we obtain 
an approximate expression for \z\ — y/ad: 



4 L J 
+ i[^^cos(2^v^)] 2 

- \[£Alsin(2gt^n)]\ (11) 

The summations in (11) involving cos(2gty/n) and sm(2gty/n) also cannot be completed 
analytically, but are the same type as those for qubit inversion in the original discussion of 
quantum revivals [10] for zero detuning, so we expect to see similar revival behavior here. 
Fig. 2 shows our analytical results for the function 2(\z\ — y/ad) in comparison with the 
complete numerical evaluation of Eq.(3) without making the X-type approximation. 

We note that both curves in Fig. 2 indicate repeated occurrences of early-stage decoherence 
(ESD). The recurring positive entanglement events are not periodic in amplitude, as found in 
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Fig. 2. The analytical and numerical results for entanglement. As expected, revivals occur, 
which arc predicted by the approximate analytical results reasonably well. Better resolution 
of the rapid oscillations is provided in Fig. 3. 



the weak-field cases discussed in [4, 5] , but are periodic in time. Fig. 3 below shows that the 
analytic approximation includes micro-ESD events that are not present in the full expression. 

Because of the exponential term in (11) the envelope of the function decreases slowly. Also, 
close to the revival regions the function makes rapid oscillations with period r = i\j(2g\fn) 
which is half the corresponding period in the inversion case. The reason is that in (11) we have 
squares of summations rather than the summations for the original inversion calculations [10]. 
As the figure shows, our analytic summations and the approximations they are based on 
work rather well, capturing all major aspects of the full numerical result. The function 
2(|z| — \fad) can suitably predict the magnitude and revival of the entanglement while not 
a perfect substitute for it. 

In conclusion, we have made a numerical and analytic examination of the entanglement 
dynamics between two qubits controlled by (or, receiving "instructions" from) quantized 
optical fields that are modelled as equal-amplitude coherent state modes. We show that pe- 
riodic ESD effects, as were predicted for very weak control fields in [4,5], are still present but 
in substantially modified form. Many very rapid appearances of the ESD effect are contained 
in the top plot of Fig. 3 and are not present in the numerical evaluation of concurrence for 
the complete two-qubit reduced p. Both the X-state simplification in format and the high- 
n approximation contribute to the differences, but the latter is much less significant when 
n^> 1, as here. A visible consequence in both plots coming just from the open-system nature 
of the coherent states is the imperfect recovery of entanglement in each successive revival 
zone. 
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Fig. 3. A more detailed plot of the results around t = 20ir/g shown in Fig. 2. Analytical 
results are for the X-type p while the numerical ones are for the original p. 



Figs. 2 and 3 have implications for further work. They show that the major features of 
the entangled qubits' response are very well captured even by the rather severe X-state 
simplification we introduced. The benefits of having analytic expressions, as provided by the 
X-state format, can be expected to be substantial in guiding and interpreting more complex 
calculations. 

In this regard, we believe it will be interesting to investigate dynamic entanglement be- 
havior under different assumptions about the state of the system. For example, the response 
to squeezed fields or fields acting at different times or fields of significantly different intensity 
or mode frequency are open for study. Similarly, the state reduction employed to reach the 
X state is not limited to two-qubit situations. These expanded topics will be the focus of a 
wider investigation. 
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